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While the coherent potential approximation (CPA) is the prevalent method for the study of 
disordered electronic systems, it fails to capture non-local correlations and Anderson localization. To 
incorporate such effects, we extend the dual fermion approach to disordered non-interacting systems 
using the replica method. Results for single- and two- particle quantities show good agreement with 
cluster extensions of the CPA; moreover, weak localization is captured. As a natural extension of 
the CPA, our method presents an alternative to the existing cluster theories. It can be used in 
various applications, including the study of disordered interacting systems, or for the description of 
non-local effects in electronic structure calculations. 
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Introduction.- Disorder, due to doping, impurities or 
structural defects, is a common feature of many materi- 
als, and may play a crucial role in determining transport 
properties. The coherent potential approximation (CPA) 
is the most common theoretical method used to study dis- 
ordered systems [l[ . While the CPA is a successful single- 
site theory of disorder it misses, by construction, non- 
local effects including Anderson localization [2[. Cluster 
extensions of the CPA, such as the Molecular CPA J3[ 
and the Dynamical Cluster Approximation (DCA) [j], 
capture non-local correlations within the cluster Q ; how- 
ever, Anderson localization is still missing. 

The Dual Fermion (DF) formalism @, originally de- 
veloped for clean interacting systems, is complementary 
to these cluster approaches. It treats local correlations 
explicitly in the "impurity" solver, and non-local correla- 
tions perturbatively. So, if a geometric series of relevant 
diagrams are included, it has the potential to capture 
localization. Here we present such a DF method for dis- 
ordered systems. 

The DF formalism is based on a set of auxiliary vari- 
ables (dual fermions) which are introduced into the path 
integral representation of the lattice partition function 
via a canonical transformation 0, Q • It maps the lattice 
onto an impurity embedded in a self-consistently deter- 
mined dual fermion lattice. The DF lattice problem is 
treated via a perturbation theory involving the DF bare 
Green function, which is the difference of the lattice and 
impurity Green functions, and the impurity full vertex as 
the effective bare DF interaction Q. These features are 
elegantly incorporated in Rubtsov's DF formalism Q. 

For systems with disorder, the DF mapping has to be 
done differently. In particular, as observable quantities 
are calculated by derivatives of the free energy, the DF 
formalism for the disorder case then needs to be con- 
structed from the disorder averaged (InZ)^^ instead of 
from the partition function Z as in Q . 



In this Letter we employ the replica method [10| to 
deal with such averaging. We extend the DF method to 
systems with disorder, and construct the DF mapping 
directly on the Green function. We demonstrate that 
our method shows remarkable agreement for the single- 
particle Green function with the results obtained from 
large cluster DCA calculations. Moreover, our disorder 
DF can account for weak localization in the conductiv- 
ity. This scheme presents a powerful alternative to the 
existing cluster extensions of CPA, with a broad venue 
of applications, including the possibility of treating both 
electron-electron interactions and disorder on equal foot- 
ing, or replacing the CPA in electronic structure calcula- 
tions. 

Method.- The simplest model of disordered electrons 
is the Anderson model. Its Hamiltonian is 

H = -J2 t^ji^Cj + h.c.) + ^ V^n„ (1) 

<ij> i 

where At = 1 sets the energy unit, and the local potential 
1^ is a site-independent random quantity, with a uniform 

1 W 

box disorder distribution, p{V) = "^®(~^ 

We focus on the disorder averaged Green function 

(G;oM),„ = -T;^(liiZ(F.,7y^O)aJw=o, (2) 

with ((•.■))q,„ = / dVp{y){...) standing for the disorder 
averaging, and "q^k is a source field. 

In the replica method, one employs the relation 

InZ = lim,„_5.o , with m being the number of 

m 

replicas. Hence, the above Green function can be writ- 
ten in terms of m powers of Z which is much more 
tractable than a logarithm. Using Grassmann func- 
tional integrals for quantum averaging, and the replica 
method for disorder averaging, one can rewrite Eq. ([2]) 
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as {Vc = Ouifca ^c"^., a the replica index 
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witli the lattice action 
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dr<(r), 
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where Wn ~ (2n+ l)7rT. Averaging over the distribution 
p{V) in Eq. (jH), we obtain 

wka i 

where W{ni) is the elastic effective interaction between 
electrons of different replicas, and may be expressed 
through cumulants < >c as Q 
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Following the DF procedure of [6[ , wc introduce an effec- 
tive single-site impurity reference problem by rewriting 
the original action as 

S = J2 5.™p[c", c"] - c»fc(A„ - Sk - v^k)Ck, (7) 

i wka 

with an effective impurity action 

S^mp = Y - M + A„,)C + W{fl,), (8) 

aw 

where A^, is a local, and yet unknown, hybridization 
function describing the interaction of the impurity with 
the effective medium. Our goal is to express the Green 
function and other properties of the original lattice via 
the quantities of such impurity problem. 

So far we have moved the local part of the lattice action 
to the effective impurity. One can go further and transfer 
the non-local part of the action to the auxiliary degrees 
of freedom, so that the original real fermions carry infor- 
mation about the local part only. In the DF scheme @ 
this is done by introducing the auxiliary dual fermions 
via Gaussian transformation of the non-local part of Eq. 
©, i.e.. 



with At 
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in which the replicated site action for site i is of the form 



^site — Simp + / ^ {c-i^fiyj + fiw'^iw 



aw 



(11) 



The inter-site coupling in Eq. ([T0| has been transferred 
to a coupling between dual fermions. This allows us to 
integrate out the real fermion degrees of freedom from 
the local site action Sl^^^ for each site i separately, i.e.. 



aw 



in which Zimp is the partition function for the replicated 
impurity system. As in the clean case formally this 
can be done up to infinite order, which makes mapping 
to DF exact. Choosing for convenience A^, = gf„ipi'^)~^ ^ 
the lowest-order of the replicated non-antisymmetrized 
DF potential [A",/f] reads as 

vTfiftJ.^] ^lY^'-'^^.^vuLfLfrn,. (13) 

where the CPA full vertex 
7"''^(u;,u;') 
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, (14) 



XoAmpiw, W') = -~g",np{w)g'^mp(w'), (15) 

with gfmpiw) = - I VcDce-^^-^^clcZ and 



xt^piw^W) = j VcDce-^-^clclclcl being the 

impurity averaged single- and two-particle Green 
functions, respectively. 

After taking the derivative with respect to the source 
field i}wki the Green function of Eq. (fTO|) reads as 



(Gfc(u;))a„ = (Au, -efe) +■ 



[(A„,-efc)5f„J«;)]^ 



(16) 



77tufe), and A^, yet to be specified. 



With such a transformation, the lattice Green function 
of Eq. (jS]) can be rewritten as 



{Gk{w)), 
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where we define the averaged DF Green function as 
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FIG. 1: (color online) First and second order DF self-energy 
diagrams. All diagrams with scattering of electrons along a 
closed loop, e.g. a and 6', vanish in the replica limit. Here, 
7 is the CPA full vertex, with conserved frequency indicated 
by dashed lines and a and P the replica. 



and = fZ, 



fwk is the non- 



interacting DF action. 

To calculate the DF Green function of Eq. (fTT]). we 
use diagrammatic perturbation theory. Here the non- 
trivial and crucial difference between the disordered and 
clean cases is that the interaction between replicas is off- 
diagonal, which puts certain constraints on the topology 
of Green function graphs. In particular, all graphs with 
closed fermion loops vanish (Fig. [Ij . This is because each 
closed fermion loop contains one free replica summation 
which gives an extra factor of m in Eq. (|17p . and thus 
equals to zero when m 
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Single-particle properties.- We first present the 
single-particle Green function for a one-dimensional 
lattice. After solving the impurity part, one 
obtains the averaged impurity Green function 

gtmp{w) = / dVp{V)- ■ --— and cor- 

responding impurity vertex 7. Then, the DF part is 
solved sclf-consistently using standard diagrammatic 
perturbation theory. Next, the real lattice Green func- 
tion from Eq. ()17|) is recalculated, and new hybridization 
function /S.{w) is constructed to parametrize the impu- 
rity problem. This is repeated until self-consistency is 
reached, namely Gd,k{w) = 0, with all local diagrams 
(e.g. diagram a in Fig. [1} being zero. 

In Fig. [5] we present results for single-particle Green 
function obtained from a fully self-consistent infinite lad- 
der diagram summation (in both particle-hole (p-h) and 
particle-particle (p-p) channels) for the DF self-energy. 
To consider the effect of non-local correlations, we com- 
pare our DF results with CPA and DCA results for clus- 
ter size Nc = 20. The local Matsubara Green func- 
tion (left panel) shows that inclusion of inter-site cor- 
relations leads to corrections to the CPA Green function. 
Both DF and DCA results show good agreement at small 
and large disorder strength. The local density of states 
(DOS) (right panel) also displays satisfactory agreement 
between DF and DCA results. Indeed, for weak disor- 



der [W = 0.25), the results from CPA, DCA and DF 
calculations are practically the same. As the disorder 
strength increases, the non-local corrections become im- 
portant (with finite momentum dependence of the self- 
energy) and the differences between the CPA and DF 
DOS are more pronounced. The DF successfully cap- 
tures such correlations by producing additional features 
which are also in good agreement with the fully con- 
verged DCA result, especially for large disorder strength 
[W = 2.0). 
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FIG. 2: (color online) The imaginary part of the local 
Matsubara Green function at T = 0.02 (left) and the to- 
tal density of states (right) for different disorder strengths: 
W = 0.25, f .25, 0.2 (4t = 1). For comparison, we present data 
obtained within CPA, a finite cluster DCA {Nc = 20) and DF 
methods. Inclusion of inter-site correlations leads to correc- 
tions to the CPA Green function (left panel) and appearance 
of additional structures at larger disorder in the total density 
of states (right panel). In each case, the DF captures the 
features of the DCA density of states and is in nearly exact 
agreement for the Green function. 

Two-particle properties.- Although the CPA provides 
a good qualitative description of single-particle quanti- 
ties, it fails to capture Anderson localization because the 
two-particle vertex docs not depend on the transfer mo- 
mentum between incoming and outgoing particles. Thus, 
the conductivity reduces to the bare p-h bubble so that 
vertex corrections are ignored 12 , 1^ . 

In our scheme, the full vertex is non-local, so we expect 
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FIG. 3: (color online) Irreducible DF p-h horizontal vertex 
F'' is calculated using parquet equation with crossing con- 
tributions from p-h "vertical" and p-p channels. The fully 
irreducible vertex is approximated by 7. F'*"'''' and F"^'^^ are 
the irreducible and full p-p vertices, respectively. 



FIG. 4: (color online) Conductivity as function of disorder 
strength using CPA and DF methods is shown. Data are 
obtained for d = 1 and d = 2 at T = 0.02. Our results show 
that vertex corrections incorporated in DF approach allow to 
capture weak localization leading to the decrease of the net 
conductivity. 



to obtain finite vertex corrections describing "weak" lo- 
calization effects IJ] . As our formalism is best converged 



on Matsubara frequency, we calculate the low tempera- 
ture dc conductivity following [l^ as 



TT 



(18) 



/3 — l/fcsT, and current-current correlation function 
Aa:x(q,T) = {jxisi,T)j^{-ci,{))). To obtain this lat- 
tice correlation function, one needs to calculate the DF 
two-particle Green function x*^ = — Xo ~ Xo-^^^Xo' '^i^h 
Xo = C^C^ [1]. For the disordered case, one has to re- 
member that in the DF vertex F'' all diagrams containing 
closed loops are zero. 

As usual, F''' is obtained from the Bcthc-Salpeter equa- 
tion F"^ =V<^ + T'^xiF'^, where F'' is the irreducible DF 
vertex in the p-h horizontal channel (c.f. Fig. [S]). To 
calculate the last quantity, we use the parquet equations 
which account for the crossing contributions from the p-p 
and the "vertical" p-h channels 0. Here the fully irre- 
ducible vertex is approximated by the impurity full ver- 
tex 7. This procedure allows us to incorporate the impor- 
tant maximally-crossed diagrams The conductivity 
can be decomposed into two parts, <j = uq -\- Act, where 
ctq is the mean-field Drude conductivity, coming from the 
bare bubble xo, and the second part Aa incorporates the 
vertex corrections. 

Our results for the dc conductivity in dimensions d ~ 1 
and d = 2 are presented in Fig. 21 The data show that 
the disorder DF method is able to capture weak localiza- 
tion with vertex corrections (vanishing in CPA) leading 
to a net decrease of conductivity. In d = 1, as the dis- 
order strength increases, the DF vertex corrections are 
more pronounced, while in d = 2 they are much weaker, 
as expected [14 1. 

Conclusions.- We generalize the DF approach to 
include disordered non-interacting systems using the 



replica method. Our formalism incorporates non-local 
disorder-induced correlations neglected in the single-site 
CPA. While in the DCA the multi-scattering effects are 
limited by the size of the cluster, our method allows us to 
treat spatial correlations on all length scales. Comparing 
our results with large-cluster DCA data, we find rather 
good agreement. This shows that the DF scheme presents 
a powerful alternative to the existing cluster theories for 
description of non-local physics in disordered systems. 
More significantly, our method incorporates finite weak 
localization corrections to the mean-field conductivity - 
the precursor effect of Anderson localization. We believe 
that the DF disorder scheme traces a clear avenue to 
study a wide variety of physical phenomena, including 
the interplay of weak localization effects and strong elec- 
tron interactions, which may be treated on equal footing 
in our method. Work on this direction and generaliza- 
tion to cluster cases [l6j are in progress. As a promising 
candidate to replace CPA, its application to study non- 
local effects in electronic structure calculations is also 
envisioned. 
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